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Abstract — Random spatial models are attractive for modeling 
heterogeneous cellular networks (HCNs) due to their realism, 
tractability, and scalability. A major limitation of such models 
to date in the context of HCNs is the neglect of network traffic 
and load: all base stations (BSs) have typically been assumed 
to always be transmitting. Small cells in particular will have a 
lighter load than macrocells, and so their contribution to the 
network interference may be significantly overstated in a fully 
loaded model. This paper incorporates a flexible notion of BS 
load by introducing a new idea of conditionally thinning the 
interference field. For a A' -tier HCN where BSs across tiers differ 
in terms of transmit power, supported data rate, deployment 
density, and now load, we derive the coverage probability for 
a typical mobile, which connects to the strongest BS signal. 
Conditioned on this connection, the interfering BSs of the i 
tier are assumed to transmit independently with probability p t , 
which models the load. Assuming - reasonably - that smaller cells 
are more lightly loaded than macrocells, the analysis shows that 
adding such access points to the network always increases the 
coverage probability. We also observe that fully loaded models 
are quite pessimistic in terms of coverage. 

Index Terms — Heterogeneous cellular networks, load-aware 
model, Poisson point process, stochastic geometry, HetNet per- 
formance analysis. 



I. Introduction 

ADVANCES in hardware and the increasing popularity of 
smartphones and tablets have led to a paradigm shift in 
the way cellular networks are accessed and consequently the 
way they are deployed. In terms of network access, focus has 
shifted from voice-oriented applications towards data-hungry 
applications such as live video streaming and symmetric video 
calls [2|. Macrocell based conventional cellular networks were 
primarily designed to provide coverage and are clearly not 
capable of accommodating this huge change in the usage 
trends [3|. One of the most promising ways to handle this 
data deluge is to increase the density of the BSs, thereby 
reducing the frequency reuse distance and hence improving 
network capacity (4). For example, a typical 3G or 4G cellular 
network already has operator-managed picocells deployed at 
the hot spots and cell edges |5 |; distributed antennas deployed 
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to eliminate coverage dead-zones (6); and low-power user- 
deployed femtocells |4); along with the existing high-power 
tower- mounted traditional macrocell BSs that guarantee uni- 
versal coverage. 

A. Related Work and Motivation 

The rapidly increasing heterogeneity and randomness in cel- 
lular networks threaten classical models, such as the hexagonal 
grid [7] or Wyner model jSJ, with obsolescence. Clearly, a 
more sensible way to model a HCN is with a random spatial 
model, where the BS locations form a realization of some 
random spatial point process j9)-fT2). Such a model captures 
the inevitable uncertainty in their locations, and tools from 



stochastic geometry 1 13 1 and point process theory 1 14 1 can be 



deployed to assist in analysis (T3) . 

For example, in an HCN, macrocells would usually follow 
a somewhat sparse point process and have high transmit 
power, whereas pico and femtocells are drawn from suc- 
cessively denser (more BSs/area) processes, and have lower 
transmit power. In such a network, a mobile user could simply 
connect to the strongest base station signal, with the rest 
of the transmitting BSs being interferers. This model was 
introduced in (16), (17) and extended in (18)-||20), and is 
surprisingly tractable: under fairly benign assumptions, the 
coverage probability could be derived in closed-form, which 
is not possible even for 1-tier networks in the hexagonal grid 
model. The model further was shown to generally agree in 
several important ways with more sophisticated industry (e.g. 
3GPP) simulations [21] and even early field deployments of 
HCNs (22). 

Despite this encouraging progress, these models lack in at 
least one important aspect, which is their neglect of network 
traffic and load. Rather, the work to date in this direction has 
assumed that all the BSs transmit concurrently all the time, 
which translates to a fully loaded (or full buffer) scenario 
resulting in pessimistic estimates of coverage and average rate. 
Although this might be justified for macrocells in peak traffic 
hours, this is not applicable for smaller cells whose smaller 
coverage areas will naturally accommodate fewer users, even 
if considerable biasing towards the small cells is introduced. 
Therefore, the main goal of this paper is to incorporate 
a notion of BS load. Those familiar with random spatial 
models will recognize that a simple independent thinning of 
the point processes will not capture the load since it may 
also turn off the serving BS, which is not allowed if the 
analysis is performed for a typical active user. On the other 
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hand, incorporating more sophisticated queueing models in the 
present multi cell scenario will render the analysis intractable 
due to the interference induced coupling in the service rates 



of various BSs |23 1, |24| . Moreover, this line of thought is not 
in the scope of the current paper since we do not focus on the 
flow level performance evaluation. The readers interested in 
flow level models can refer to fl25) , (26). With our main focus 
on the downlink coverage evaluation, we propose a middle way 
whereby we conditionally thin the interference field predicated 
on a connection to a typical active user, and we are able to 
maintain acceptable tractability with a realistic model of BS 
loading. 

B. Contributions and Outcomes 

The main contributions of this paper are as follows: 

1) Tractable Load Model for K-Tier HCNs: In Section \U\ 
we incorporate a notion of BS load in a general K-tier random 
spatial model for HCNs. For an HCN where BSs across tiers 
differ in terms of their transmit power, supported data rate 
and deployment density, we assume that a typical mobile 
connects to the strongest BS in terms of received power and 
conditioned on this connection, the i th tier interfering BSs 
transmit independently with a probability p. L , which models 
the load. These BS activity factors {p,^ may vary significantly 
across the tiers due to different coverage areas of each tier. 

2) Coverage Probability: We derive exact expressions for 
the coverage probability of a typical mobile user in both open 
and closed access HCNs. Since these expressions involve an 
infinite summation, we also derive a set of upper and lower 
bounds that can be made arbitrarily tight with a finite number 
of terms. These bounds also give insights into the number 
of terms of the infinite summation required to approximate 
the coverage probability such that the approximation error is 
within some predefined limit. 

3) Design Insights: This paper provides some potentially 
useful design insights for HCNs. First, we study the effect 
of proposed "conditional thinning" on the coverage footprints 
of various tiers and show that this effect can be understood 
in two equivalent ways: i) thinning of interference, and ii) 
biasing of the typical mobile towards its serving BS. While 
the former is a direct result of thinning, the latter is an indirect 
consequence of the expansion of the coverage regions in the 
thinned interference field. 

Second, our analysis sheds light into the effect of adding 
new tiers to already existing HCNs. In particular, we derive 
an exact condition under which the addition of a new tier 
to a general i^-tier HCN will increase the overall coverage 
probability. A relevant special case is the addition of small 
cells to existing macrocell networks, where we show that 
in the interference limited regime the overall open access 
coverage probability increases if the load on small cells is 
smaller than that of macrocells, which is a typical operating 
scenario because of the smaller loads handled by small cells. 
This is a strong rebuttal to the viewpoint that unplanned 
infrastructure might bring down a cellular network due to 
increased interference. 

Third, we show that the coverage probability for a general 
K-tier interference-limited open-access network is invariant 



to changes in the power and deployment density when all the 
classes of BSs have same loads and target SINRs. Furthermore, 
this coverage probability is also the same as that of a single tier 
network with the same target SINR and the same BS activity 
factor. 

II. System Model 

We model a downlink heterogeneous cellular network with 
K classes (or tiers) of BSs. For notational simplicity, we 
denote the set {1, 2, . . . K} by K,. BSs of the i th class transmit 
with power Pi, have a target SINR of fa and are assumed to 
form a realization of an independent homogeneous Poisson 
Point Process (PPP) $i with density A^. Such a model seems 
sensible for user deployed BSs such as femtocells but is 
dubious for the centrally planned tiers such as macrocells. 
Nevertheless, the difference is not as large as expected and 
PPP assumption for macrocells is shown to be about as 
accurate as the grid model when compared to an actual 4G 
network in [ flT) . More recently, |27) has validated the PPP 
assumption for certain cities using tools from spatial statistics. 
Furthermore, this model is likely even more sensible for K- 
tier HCNs due to the increased uncertainty in the deployment 
of lower tiers (smaller cells). We will comment more on the 
accuracy of this assumption in the light of the proposed load 
model in the Numerical Results Section. 

Without loss of generality, we perform analysis on a typical 
mobile user located at origin, which is made possible by 
Slivnyak's Theorem fl3) . For cell association, we consider 
the max-SINR connectivity model, where a mobile user con- 
nects to the BS that provides highest downlink SINR. It 
should be noted that this model is the same as the max- 
power connectivity model where a mobile connects to the 
BS that provides highest downlink power. Since HCNs are 
typically interference-limited (28), we ignore thermal noise 
for notational simplicity. However, as would be evident from 
the analysis, this assumption can be relaxed without much 
extra work. To model the wireless channel, we consider a 
standard distance based path loss with exponent a > 2 along 
with Rayleigh fading. Hence the received power at a typical 
mobile from a BS located at point x E $i can be expressed as 



Pih, 



where h x ~ exp(l) and 



is the distance 



based path loss. General fading distributions, e.g., log-normal 
shadowing, can be incorporated using techniques developed 



in 1 29 1 at the cost of tractability. Assuming to be the set 
of kF^ tier interfering BSs (possibly thinned version of 
the downlink SIR at the typical mobile user when it connects 
to the BS located at point y G $; is: 



SIR(y) = 



PihvWvt 



T,£=iY, xe Zu P khx\\x 



(1) 



A. Modeling Base-Station Load 

In this A"-tier random spatial model, we now incorporate 
network "load" perceived by each BS as the likelihood of its 
transmission at a randomly chosen time instant. This can also 
be visualized as the BS activity factor, formally defined as the 
fraction of time for which a BS transmits. 
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Relationship of BS activity factor with number of active 
users: A BS is inactive in a particular resource block, e.g., 
time-frequency resource block in LTE [30|, if there is no active 
user scheduled. This can be due to an over provisioned system 
or a momentary lull in traffic due to the bursty nature of data 
access. Clearly, this model characterizes the load on each BS 
in terms of the total number of active users served by that 
BS at a random time instant. In the context of Orthogonal 
Frequency Division Multiple Access (OFDMA) if a particular 
BS experiences high load, it will utilize more frequency time 
resources and hence the probability that a user is scheduled 
in a particular frequency time block increases. Therefore, the 
load perceived by a BS is directly related to the likelihood that 
an arbitrary resource block is utilized and hence is related to 
the BS activity in that particular block. 

Temporal and spatial correlation in BS activity factors: 
In general, there is both temporal and spatial correlation in 
the activity factors of different BSs. Temporal correlation is 
induced across neighboring BSs by the mobility of users, i.e., 
if a user is associated to a particular BS, the likelihood of 
neighboring BSs transmitting at a future time instant is slightly 
higher. Spatial correlation is induced by interference and 



and the total received power at the origin from the set of active 
BSs as: 



traffic/load patterns ]23[ , 1 24 1. To understand this, consider two 
neighboring BSs. When the first BS transmits, it increases net 
interference experienced by the second BS and hence reduces 
its data rate. As a result, the second BS now takes longer to 
transmit same amount of data than it would have taken if the 
first BS was not transmitting. Therefore, the activity factors of 
these two BSs are positively correlated. However, modeling 
the exact nature of these correlations is beyond the scope 
of the current paper and we assume the BS activity factors 
to be independent. Although the spatio-temporal correlations 
haven't yet been modeled for this exact problem, it is worth 
noting that they have been handled in some related setups, 
e.g., the effect of spatio-temporal correlations of interference 



on coverage is discussed in [31], [32 1. 



B. Proposed Load Model and Mathematical Preliminaries 

We assume that a typical mobile connects to the strongest 
BS in terms of received power and conditioned on this 
connection, the interferer belonging to the i th tier transmits 
independently with a probability pi and is idle with a prob- 
ability 1 — pi. This conditioning makes it harder to analyze 
this system model since we do not have a priori knowledge 
about the serving BS and hence it is not possible to isolate the 
interference field. To overcome this, we partition each tier <3> m 
independently into two sets of BSs ^f m and A m , where ^ m 
and A m are both independent PPPs with densities p m X m and 
(1— Pm)A m . The set ^ TO represents the set of active BSs of tier 
m with the possibility of one of them being a serving BS, and 
A m represents the set of idle BSs of tier m with an exception 
that it could also contain the serving BS since partitioning was 
done independently. The advantage of this partitioning is that 
the interferers are confined to the set ^> = UmGK:^' m - ^ or 
ease of notation, we define the maximum signal strength from 
a set of nodes A as 



£ £ PM\x\ 

i=l 



(3) 



which denotes the net interference power if ^ does not include 
the serving BS and the interference plus signal power if it 
includes the serving BS. From the definition of M(\Pi) and 

/, it is easy to see that 1 ^ i^m(-1) < @y = ^ on ^y ^ no 
active BS in the set itj can connect to the mobile. Similarly, 
1 ^Mlpl < fi\ = i on iy if no BS in the set A; is able to 
connect to the mobile. The second event is defined to cover 
the possibility that a serving BS may lie in the set Aj. Using 
these two events, we will now define the coverage probability 
of a typical mobile at the origin. We note that a mobile will 
be in outage (not in coverage) if none of the BSs in the whole 
network provides SIR that is greater than the corresponding 
target for that tier. 

Definition 1 (Coverage Probability). Coverage probability, P c , 
can be formally defined as: 



1 - E 



- M(#) 



< A 



i(^<A 



(4) 



M{A) = sap P A h x \\x\\ 



(2) 



For this definition, we implicitly assumed an open access 
network where a mobile user is allowed to connect to any 
BS in the network without any restrictions. Another possible 
access strategy is closed access or closed subscriber group 
strategy in which a mobile is allowed to connect to only a 
subset B C JC of all the tiers. Coverage probability for closed 
access is also given by |4} with the only difference that the 
product is over the set B instead of K,. 

For tractability, we assume that the target SIR thresholds 
Pi are greater than dB, i.e., (3i > 1, V i. This is in fact 
the case for a large fraction of mobile users and only a few 
edge users might violate this assumption. Moreover, in the 
Numerical Results Section we show that the results derived 
under this weaker assumption hold down until around — 2dB 
which covers a large fraction of cell edge users as well. This 
assumption has also been validated earlier for the fully loaded 
K-tiei HCN in [ 17]. The reason why this assumption is helpful 
is because it ensures that at most one BS in the active set \& 
meets the target SIR requirements for a typical mobile user. 
Refer to fT7) for a detailed discussion on this assumption and 
its application in coverage analysis of a fully loaded X-tier 
HCN. 

C. Coverage Regions 

Before going into the detailed analysis of coverage prob- 
ability, it will be useful to understand the effect of the 
proposed load model on the coverage footprints of various 
BSs. Consider a realization of a three tier HCN in Fig. [T] We 
first plot the coverage regions assuming a fully loaded network 
by tessellating the space according to max-SIR connectivity 
model in the left figure. Clearly, this plot does not resemble 
a classical Voronoi tessellation due to the differences in the 
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transmit powers of BSs across tiers. Moreover, it should be 
noted that the "cell edges" are not as sharp in reality due 
to fading and shadowing, which are averaged out for these 
illustrative plots. The effect of incorporating the proposed 
load model on coverage footprints can now be understood 
in two equivalent ways: i) thinning of the interference field 
conditional on the connection of a typical mobile to its serving 
BS, where the original coverage regions corresponding to the 
inactive BSs are removed to highlight conditional thinning 
(middle figure), ii) biasing of a typical mobile towards its 
serving BS relative to the new cell edge defined by the set 
of active BSs (right figure). While the former is a direct result 
of conditional thinning, the latter is an indirect consequence of 
the expansion of coverage regions in the thinned interference 
field. 

III. Coverage Probability 

This is the main technical section of this paper where we 
derive the probability that a typical mobile is in coverage under 
the system model introduced in the last section. We first derive 
coverage probability for an open access network, from which 
the results for closed-access immediately follow. 

A. Exact Expression for Coverage Probability 

We start by stating the Laplace transform of /, i.e., £i(s) — 
E [cxp(— si)], in Lemma [TJ which will be useful in the 
derivation of coverage probability. The proof is given in flT) . 

Lemma 1. The Laplace transform of I can be expressed as: 
CM = exp (s^C(a)J2Pi^ p A , (5) 



where C(a) is given by: 

C(a) = 



27r 2 csc 



(6) 



The following Lemma deals with fractional moments of in- 
terference and is the main technical result required to evaluate 
the coverage probability for this model. 



Lemma 2. Let ^i denote the set of active transmitters of tier 
i and Si — Pi/ (1 + fti). Let I denote the total received power 
from the BSs in the set ^ and for notational simplicity define 

T = 1 f max < i) r 2 /<*. Then 



\ ie/c 



E [T m ] 



m\g(m) 
(-A) m % 



where 
g(m) 



r(i 



B 7rT(l+j) 
rj y(1 + ( m+1 ) 2 

V ' a 



and 



A = nr(l + -)J2(l-Pi)^iPi 



i > j i ' 



ie/c 



K 



T 1 = C{a)Y J PihPC 



(7) 

(8) 

(9) 
(10) 



z=i 



The hypergeometric function is denoted by 2Fi{a,b,c 1 z) 
r(c) ri t b - 1 (i-t) c -*>- 1 
Jo 



-dt 



r(6)r( c -b) Jo (i-tz)» 

Proof: See Appendix |A] ■ 
Using these Lemmas, we now derive the main coverage 
probability result. 

Theorem 1 (Open Access). The downlink coverage probabil- 
ity for a typical mobile user in a K-tier open access network 
assuming fii > 1, V i, is 



E PiKPi p, 

7T IG/C 



2/ao-Va 



EIip^p, 



,2/a 



Proof: The coverage probability is given by Q. Since 
the point processes Aj and the corresponding fading random 
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variables are independent, conditioning on the common in- 
terference, we can move the expectation inside the product. 
Hence 

A 



1-P B 



E 



mi 



M(*) 



< ft E [1 (M(A 4 ) < ft/)] 



(12) 

where the inner expectation is with respect to the inactive 
transmitter sets. We first simplify this inner expectation as 
follows: 



E [1 (M (A() < ft/)] 



E 



( = } E 



(b) 



Y[ i (Pih\\x\\- a < ft/) 
n (l-expt-ft/^/iMr)) 

•eA, 

cxp ( -(1 - Pi )Xi f exp (-ftfT 1 / 



dx 



exp 



-(l-ft)Aift a r«PfnT 1 



(13) 
(14) 

) 

(15) 
(16) 



where (a) follows form the fact that fading is Rayleigh 
distributed, i.e., h ~ exp(l), (fe) follows from the probability 
generating functional (PGFL) of PPP fl3[ and (c) follows 
from some algebraic manipulations to reduce the integral to a 
Gamma function. Now recalling the expression of A given by 
(|8), we can write: 

Af(*i) 



1 — P c = E 



1 max ■ 

I ie/C 



< / cxp(-A/- 2 /") 



(17) 



Using the Taylor series expansion of exp(— x), exchanging the 
infinite summation and expectation^] 



ml 



m=0 



\ i£K Si 

The summation can be split as: 



iL^</)r^ 



1 — P c = r [ max ■ 

V ie/c 



< / 



H [ r ' n ] ' (18) 

m— 1 



/?? 



The term 1 — P ^maxi < jj [ s me coverage probability 

in a fully loaded heterogeneous network where the m-th tier 
density is p m A m . This is derived in fl7) and is given by: 



A' 



l-p(maxffi)</! = 



i=l 



2/a R -2/a 



2/c 



(19) 

Using Lemma [2] to evaluate E [T m ], we obtain the result. ■ 
We note that the expression of coverage probability involves 
infinite summation over the sequence g(m). Therefore, we first 
show that the infinite summation converges by showing that 
|g(m)| — > as m —> oo. Observe that: 



\9(m)\ < [ - 
V 



< 



(A/vY 



'The average of the series is absolutely convergent. 



{Ah) 



r 2m ~i 
I a I 



0, (20) 



where the limiting argument follows from the fact that the 
sequence of the form x n /n\ — > 0. In addition to proving that 
the series converges, this upper bound on |ff(m)| also sheds 
light on the behavior of the sequence g(m). If A/rj < 1, the 
bound decreases monotonically with m and hence it is enough 
to consider only a few significant terms to closely approximate 
the infinite sum. However, if A/rj > 1, especially if A/rj 1, 

the upper bound first increases until \— ] < {A/rj) r^ar 1 and 
decreases thereafter. Therefore, the number of significant terms 
of g{m) required to approximate the infinite sum would be 
higher. It can be easily shown that A/rj < 1 for all choices of 
system parameters when the activity factor of each tier satisfies 
the following condition: 

1 



Pi > 



1 



C(a)fi /a [nT(l + 2/a)p 

For ft = 1 and a = 4, this value of pi comes out to be w 0.36. 
Therefore, the infinite sum can be tightly approximated by 
the first few significant terms of g{m) in most operating 
scenarios. We will comment more on the convergence of g{m) 
and the number of terms required to tightly approximate the 
coverage probability later in this section and in the Numerical 
Results Section. We now provide the exact expression for 
the coverage probability in a closed access network in the 
following Theorem. We recall that that coverage probability 
in closed-access is given by Q with the only change that 
the product is over B instead of JC. By definition, coverage 
probability in closed access is less than that of open access. 
Using this definition, the proof proceeds exactly same as that 
of Theorem [T] and hence is not provided. 

Theorem 2 (Closed Access). The downlink coverage proba- 
bility of a typical mobile in a K-tier closed access network 
where a mobile is allowed to connect to B C /C tiers assuming 

ft > 1, V i, is 



(21) 



Pc = 



C{a) 



^r K r, \ P 2/c 



- 9c{m), 



(22) 



where g c {m) and the corresponding expression for A are given 
by Q and ([8]), respectively, with the only difference that the 
summations defined over set K, are now over set B. 

We conclude this discussion with a note that the proof 
technique introduced in this section is quite general and can be 
used to study variants of the load model introduced in the last 
section. For example, if the network is modeled such that it has 
a predefined set of BSs that are active and a typical mobile is 
allowed to connect only to the inactive set, it is easy to observe 
that the coverage probability under open access assumption 



is given by P c = 1 — E 



n; 



A 



M(Aj) 



<ft 



From the 

proof of Theorem[l] this corresponds to 1 — E[cxp(— AI~ 2 / a )] 
and can easily be evaluated following the proof technique of 
Theorem [T] The same argument can be extended to the closed 
access case as well. 
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B. Special Cases of Interest 

We now use the results derived in this section to study 
some special cases and compare the system performance with 
already known results for fully loaded system. First, we note 
that for a fully loaded system, the value of A — and hence 
g(m) = g c (m) = 0, V m. Therefore, the coverage probability 
in this case can be expressed as the following Corollary of 
Theorems Q] and [2] 

Corollary 1 (Fully Loaded). For a fully loaded system, i.e., 
Pi = 1 V i, the coverage probability in open access is given 
by: 

E KPf' a p; 2/a 



7T i£K 



1/a 



(23) 



which is the same as Corollary 1 in /[IT} . The coverage 
probability in closed access is also given by ( I23J I with the 
only difference that the summation over the set IC is now over 
set B. 

For a single tier open access network, the coverage proba- 
bility derived in Theorem [T] can be simplified and is expressed 
as the following Corollary. 

Corollary 2 (Single Tier). The coverage probability for the 
single tier open access network with BS activity factor p is 



P, = 



7r/3- 2 /° 
C{a) 



9(m), 



(24) 



where the terms — and — appearing in the expression of g(m) 



given 



by Q, 



.4 

n 

B 
V 



7rr(l + f)(l-p) 



C(a) P pi 



1 



(m+l)2 



1 



(25) 



(26) 



C(a)0S(l + /3)^ 

Remark 1 (Scale invariance of a single tier network). From 
Corollary [2] we note that for any BS activity factor p, the 
coverage probability in a single tier open access network is 
independent of the BS density A and transmit power P. This 
is henceforth referred to as "scale-invariance" of cellular 
networks to changes in the BS density and their transmit 
powers. 

Remark [T] is a generalization of a similar result derived for 
fully loaded networks in p7[ , which can easily be seen from 
Corollary [T] In addition to single tier networks, it was also 
observed in flT) that the general fully loaded open access 
multi tier networks also exhibit scale invariance if the target 
SIRs for all the tiers are the same. This can also be easily 
deduced from Corollary [T] Motivated by this observation, we 
study the coverage probability for our proposed load model in 
open-access multi tier networks under the assumption that the 
target SIR is the same for all tiers in the next Corollary. 

Corollary 3 (Coverage Probability: if -Tier with same j3). The 
coverage probability for a K-tier open access network under 
the proposed load model assuming target SIRs to be the same 



(= (3) for all the tiers is given by ( j24| >, with the difference that 
the term ~ appearing in the expression of g(m) given by Q 
is: 



V 



(27) 



and — appearing in |7]) is given by ( |26| l. 

Remark 2 (Scale invariance of K-tier HCNs with same j3). 
From Corollary^ we note that the coverage probability for K- 
tier HCNs is not scale invariant in general, even when target 
SIRs of all the tiers are the same. However, the invariance 
property does hold when the BS activity factors of all the 
tiers are the same. Interestingly, the coverage probability in 
this case is same as that of a single tier network given by 
Corollary [2] 

To understand this remark, we consider the following simple 
example. 

Example 1 (Scale invariance in a 2-tier HCN). Consider a two 
tier network with BS activity factors p\ and P2- If Pi < Pi, 
increasing the density of the first tier leads to a higher increase 
in the intended power due to the higher likelihood of having 
a closer tier-1 BS as the serving BS but a relatively smaller 
increase in the interference power. The coverage probability 
in this case is expected to increase. On the other hand, if 
Pi > Pi, increasing the density of tier-l BSs leads to higher 
increase in the interference power as compared to the intended 
power, leading to a decrease in the coverage probability. The 
two effects cancel each other when the activity factors of the 
two tiers are the same. 

We now extend this result and derive exact condition under 
which the addition of (K + l) th tier won't affect (or will 
improve) the coverage of the existing if -tier network. We 
again assume same target SIR for all the tiers. The result is 
given in the following Corollary. 

Corollary 4 (Same j3: Effect of adding (K + l) th tier). The 
overall coverage probability increases with the addition of the 
(K + l) th tier if the load on the new tier satisfies: 



K 

PK+1 <^Pl 



1/a 



2/a' 



(28) 



decreases if the inequality is reversed and remains the same 
if d28| holds with equality. 



Proof: From Corollary [3] we note that the only term in 
the coverage probability expression that will change with the 
addition of a new tier is A/77. ^ can t> e expressed as: 



A = TrtXl + f ) 
7] C(a)0 a 



1/a 



K 

.1=1 l^i=l A % r i 



1/a 



Defining effective load on a if- tier network as: 



(K) 



K 



Pi-. 



1/a 



1=1 Ej=l ^i-Pj 



1/a ' 



(29) 



(30) 
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A/r) can be expressed as: 
A vrr(l 



{K) 

PeS 



?7 C(a)f3 2 / a 



(K) 
PeS 



(31) 



which is the same as ( |2"5j ) for the single tier coverage result 
derived in Corollary [2] From this equivalence, it follows that 
the coverage probability is a decreasing function of p e g. 
Therefore, if the addition of the new tier leads to lower 
effective load on the network, the coverage will increase. This 
can be shown to be the case when d28ll holds as follows: 



PeS 



E 

i=i 



X,P, 



2/a 



Pi 



E 



K+l , p 2/c 
1 A i^i 



(32) 



< 



Y,?=iPih p i 



Hi=l YlnLl 



(K) 
■PeS ■ 



(33) 

(34) 



The other two results follow using the same argument. 



C. Bounds on the Coverage Probability 

Evaluation of the exact expression of the coverage probabil- 
ity requires an infinite summation. Although we have argued 
that the summation can be tightly approximated by considering 
only a first few terms, we haven't yet provided a formal 
method to determine the exact number of terms required such 
that the approximation error is within predefined limit, say e. 
Interestingly, this can be achieved as a by-product of the set of 
bounds we derive in this section that can be made arbitrarily 
tight. The idea is to use the following identity of exp(— x). 

Lemma 3. For x > and to > 0, 

2m— 1 / \ a 2m , x ,- 

(-«) / , _n / V- (-aO 



E 



< exp(-x) < ^2 



(35) 

i'. 

i=0 i=0 

Proof: The proof follows from induction. Since 
exp(— 0) = 1 and EST* ^~7r~ = 1' ^ suffices to prove 

that 4- E?™o _1 ^tt exp(— x) < 0, which follows from 

the upper bound when m = rn — 1. ■ 
Using this identity in the proof of Theorem [T] results in the 
following bounds. 

Lemma 4 (Bounds on Coverage Probability). For m > 0, 

the coverage probability for the proposed load model can be 
bounded as 

E X \ t-)2/q q — 2/a 
PiXiPi' Pi 2m- 1 

>(0<Pc- 



2m 

E: 

i=i 



2/a 



< - E -9» 

i=l 

(36) 



Clearly, these bounds can be made arbitrarily tight by in- 
creasing the value of m. Interestingly, these bounds are closely 
related to the exact expression of coverage probability derived 
in Theorem [T] In particular, the upper and lower bounds are 
derived by truncating the infinite sum over g(m) at odd and 
even number of terms, respectively. Therefore, these bounds 
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Fig. 2. Plot showing the convergence of the series J\ g(i) for various BS 
activity factors in a single tier network with 0=1, 



provide a direct way to find the number of terms of g(m) 
required to ensure an approximation error within a predefined 
limit e, which is equal to M e , where M e = min |.g(m)| < e. 

m 

We will use this observation in the study of the convergence 
of infinite sum over g(m) in the Numerical Results Section. 

We conclude this section by noting that some terms of the 
sequence g(m) can be expressed in closed form, leading to 
closed form bounds for the special case when a = 4 and 
m = 2. The bounds in this case depend only on the first two 
terms of g (to) that can be expressed as: 

-A \2 V^Ap^^j (3?) 




i - — E x iPi p iiPi 5 - csc-^yr+A)) 

^ i=l 



IV. Numerical Results 



(38) 



Since most of the analytical results derived in this paper are 
fairly self-explanatory and do not require separate numerical 
treatment, we will provide only those results which help in 
validating key modeling assumptions or help better visualize 
certain important trends. 

A. Convergence of Infinite Sum 

We study the convergence of the infinite sum appearing in 
the coverage probability expression in Figs. [2] and [3] Fig. [2] 
plots the truncated series E2=i as tne function of to for a 
single tier network and hence gives insights about the number 
of terms required until the series converges. To understand the 
trends, recall that the ratio A/rj decreases monotonically with 
the activity factor p. Therefore, the number of terms required 
for the series to converge are higher when the BS activity 
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Transmission Probability 



Fig. 3. Number of terms of the sequence g(m) required as the function of 
the transmission probability of the lowest tier for e = 10 — 8 . 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
Transmission Probability (p) 



Fig. 4. Coverage probability as a function of transmission probability is a 
single tier network (/3 = 1 and a = 3.8). 

factor is lower. Moreover, for the case when A/rj > 1, i.e., 
p = .25, the series first increases until a certain point and 
then decreases and finally converges to its limiting value. This 
trend has been discussed in detail earlier in the paper when we 
proved the convergence of the infinite sum. To provide an idea 
of the number of terms required such that the approximation 
is within e of the exact value, we plot the number of terms M e 
for various scenarios in Fig. [3] We again note that the number 
of terms required are reasonably small unless the transmission 
probability of some tier is extremely small. 

B. System Model Validation 

1 ) Comparison with the fully loaded system: After gaining 
insights into the behavior of the coverage probability expres- 
sion, we now use it to highlight the importance of the proposed 
model by comparing the coverage results of a single tier 
network with those of a fully loaded system in Fig. [4] Although 




Target SIRP (in dB) 



Fig. 5. Comparison of the coverage probability of the PPP and grid models 
with the actual BS locations of macrocells. The second tier is modeled as 
PPP in all three cases (K = 2, A 2 = 2Ai, P = [1,0.01], L = 40 X 40 
Km 2 , a = 3.8, p = [0.8, 0.6], @i = fa = /?)■ 



a huge difference in the coverage guarantees was expected for 
very low BS activity factors, it is indeed interesting that the 
coverage estimates assuming full load are quite pessimistic 
even for reasonably high load scenarios, such as p — .7 — .8. 

2) Comparison with an actual AG deployment: After high- 
lighting the importance of the proposed load model, we now 
validate the PPP model used for the BS locations in this paper. 
While this model seems sensible for the small cells, especially 
the ones driven by unplanned user deployments, such as 
femtocells, it is dubious for the centrally planned tiers such as 
macrocells. Therefore, with a special focus on the macrocells, 
we consider three location models for a two tier HCN: i) 
macrocell locations modeled as a realization of a PPP, ii) 
macrocell locations modeled as a hexagonal grid, iii) macrocell 
locations drawn from an actual 4G deployment over 40 x 40 
Km 2 area fll) , (lT\ . The second tier is modeled as a PPP in 
all three cases. The numerically evaluated coverage probability 
results for all these models along with the analytical results 
of the proposed load model and the fully-loaded PPP model 
are presented in Fig. [5] We first note that the proposed PPP 
model is about as accurate as the grid model when compared 
to the actual 4G deployment, with the grid model providing 
an upper bound and the PPP model providing a lower bound 
to the actual coverage probability. This is consistent with 
the conclusions of pT[ , p7| , which focus on fully-loaded 
cellular models in single tier and multi tier cellular networks, 
respectively. Second, we note that the analytical results derived 
for the proposed load model are accurate down to about — 2dB 
even though they were derived under the assumption that the 
target-SIR is greater than OdB for all the tiers. Since this 
covers most of the cell edge users as well, the proposed 
analytical results are reasonably accurate in the operational 
regime of the current cellular networks. Third, we note that 
the fully-loaded model provides a very loose lower bound to 
the actual coverage probability, thereby highlighting again the 
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Fig. 6. Comparison of the derived theoretical results with detailed system 
simulation accounting for actual load factors resulting from actual coverage 
regions (K = 2, P = [1,0.1], ft = j8 2 = OdB, Ai = A 2 , Mfl B = 20, 
a = 3.8). 



importance of the proposed load model. 

3) Comparison with a detailed simulation: The following 
two assumptions were made to facilitate analysis: i) the activity 
factors are the same for all the BSs of a particular tier, 
and ii) the activity of each BS is independent of the other 
BSs. We now validate these assumptions by comparing the 
analytical results with a detailed system simulation. For this 
comparison, we consider a simulation setup consisting of a two 
tier HCN, with the BSs of each tier modeled by independent 
PPPs. The user locations are also modeled by an independent 
PPP with density A u . As in the proposed model, each user 
is associated with the BS that provides the best received 
signal strength. From this, we calculate the actual load being 
served by each BS in terms of the number of users, which we 
denote by N x . for a BS located at a;,. Assuming the number 
of orthogonal resource blocks, e.g., time-frequency resource 
blocks in LTE pO] , to be M, the activity factor of a BS in 
each resource block can be expressed as p Xi — N X JM as 
discussed in Section [D] To keep the setup simple, we consider 
the regime where the probability of having N x . > M for 
any BS is small and whenever it happens, the activity factor 
for that BS is assumed to be 1. For this setup, the coverage 
probability results are presented in Fig. [6] 

For a meaningful comparison of this simulation result 
with the analytical results, we first need to find analytical 
expressions of the activity factors pi as a function of the user 
density X u . For this, we leverage Corollary 2 of p7) , where it 
is shown that the fraction of users served by j th tier is given 
by: 



Fig. 7. Coverage probability in a two tier network as a function of A2 
03 = [1, 1], P = [1, .01], Ai = 1, pi = .6 and a = 3.8). 



(39) 



Using this result, the average number of users served by 
a j th tier BS (average load) is &Nj. Therefore, assuming 
M resource blocks, the activity factor in a randomly chosen 



resource block is: 



Pi 



1 X U N, X u {Pj/Pj)' 



M X, 



M 



(40) 



We use this analytical result for the the load factors in the 
coverage probability results derived in the paper to plot the 
analytical results as a function of the user density in Fig. [6] 
Comparing this result with the numerical result obtained form 
a detailed simulation, we note that the two are reasonably 
close, especially relative to the previously known results for 
the fully-loaded system. This validates the two assumptions 
mentioned in the starting of this discussion. 

C. Scale Invariance and Effect of Adding Small Cells 

We now consider a two tier system and plot the coverage 
probability as a function of the density of second tier for 
various BS activity factors in Fig. [7] The target SIR is fixed to 
be the same for both the tiers. We first note that the network is 
invariant to the changes in density when pi = p 2 as discussed 
in the last section. More importantly, we note that the coverage 
probability increases with A2 when the second tier BSs are 
less active than the first tier. This is an important result from 
the perspective of small cells, which are generally less active 
than macrocell BSs. Therefore, the coverage probability of 
the network should increase with the addition of small cells 
in this regime. This is a strong rebuttal to the viewpoint that 
unplanned infrastructure might bring down a cellular network 
due to increased interference. On the other hand, if a tier 
of BSs is added which is more active than the macrocells, 
the coverage would decrease, although this case seems pretty 
unlikely given the high load handled by the macrocells. 

D. Open vs Closed Access 

So far in this section we have only studied open access 
networks, where a mobile user can access any BS in the 
network. We now study the effect of closed access on the 
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Fraction of second tier BSs in open access 
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Fig. 8. Coverage probability in a two tier network as a function of the 
fraction of the second tier BSs in open access (/3 = [1, 1], P = [1, .01], 



A 



(<0 



lOAi, pi = 1 and a = 3.8). The density of second tier BSs in open 



access A 



(o) 



i-/ 



where / is the fraction of BSs in open access. 



coverage probability, with a focus on a particular scenario of 
interest where a certain fraction of BSs of a particular tier 
are in closed access while the others are in open access. This 
scenario is especially important in the current HCNs, where 
the access permissions of a particular small cell might be 
different for different set of users. It is easy to argue that 
this scenario can be visualized as a special case of the general 
model developed in this paper. To understand this, assume that 
a fraction / of BSs of a certain tier are in open access - we 
assume that a BS is in open or closed access independent of the 
other BSs. Therefore, the density of BSs in open access Aj 
and closed access A^ can be evaluated from the following 
two equations: 



/ = 



A 



(o) 



Af o) + tf> 



A^A^ + A^. 



(41) 



(42) 



Now this tier can be divided into two tiers, one with density 
Xf , which is in open access, and other with density A^ , 
which is in closed access - both form independent PPPs. 

To study this scenario in detail, we consider a two tier HCN, 
where the first tier is in open access and fraction 1 — / of BSs 
of the second tier is in closed access. For this scenario, the 
coverage probability as a function of / is presented in Fig. [S] 
for various load scenarios. The results confirm the intuition 
that the gap in open and closed access results reduces when the 
value of / is increased. More interestingly, the gap is smaller 
when the second tier BSs are lightly loaded. This implies that 
the effect of interference due to closed access small cells on 
coverage probability is negligible if there are enough small 
cells in open access. 

V. Conclusions 

In this paper, we have incorporated a flexible notion of BS 
load in random spatial models for K-tier HCNs by introducing 



a new idea of conditionally thinning the interference field, 
conditional on the connection of a typical mobile to its 
serving BS. We have shown that this conditional thinning 
is a natural way of modeling different levels of load on 
different types of BSs arising mainly from the differences in 
their coverage footprints. We observe that the fully loaded 
models are extremely pessimistic in terms of coverage, and the 
analysis shows that adding lightly loaded access points (e.g. 
pico or femtocells) to the macrocell network always increases 
coverage probability. 

This work has numerous extensions. Firstly, it is important 
to exactly model the temporal and spatial correlation in the 
BS activity factors and develop tools to incorporate it in the 
framework developed in this paper. Secondly, it is important 
to understand the coverage and rate trends for this conditional- 
thinning using more realistic spatial models that include inter- 
point interactions, such as Gibbs models ETJ. Another related 
idea is to define BS load from queuing perspective and then 
study spatial and temporal dynamics simultaneously. Since this 
work focused only on the downlink, an interesting extension 
is to develop a similar framework for uplink possibly using 
some of the ideas recently developed in J33J. 



Appendix A 
Proof of Lemma[2} 

Being consistent with the definition of T, we note that: 



T m = 1 max ■ 



(43) 



To proceed with the proof, we represent / 2m /" j n terms of 
T(x) as: 



J — 2m/ a 



1 



T(2m/a) J 



e- sI s- 1+£ 



"ds, to > 1, (44) 



where T(x) is the standard gamma function. Using this rep- 
resentation of J- 2m / a we can express E[T m ] as: 



E 



e- sI s- 1+± 



(45) 

Using Fubini's theorem, we can exchange the expectation and 
the inner integral to obtain 



1 



"E 



; --lLax^l</ 
V ietc Si 



(46) 

Under the assumption > 1, V i, we know that only one 
BS in the whole network can establish a downlink connection 
with a typical mobile. Hence, 



1 max 



>/U^^l(SIR(x)>ft), (47) 



i=l 



where SIR(x) is the received SIR when a typical mobile is 
camped to the BS located at x 6 ^i. Using this expression, 
the expectation term of (|46) can be written as: 

E [ e -i( mf S</ 
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K 



e- sI l(SIR(x) > Pi 



(48) 



From Lemma [T] we know the Laplace transform of total 
interference and hence the first term in the above expression 
can be directly written as: 

E [er sI ] = exp (^s 2 / a C(a) £>Ajif /o, J . (49) 



To evaluate the expectation in the second term of ( |48| i, we first 
denote the effective interference as I' = I — Pih x \\x\\~ a and 
note that the Laplace transforms of I and I' are the same. The 
expectation can now be simplified as: 



E 



e' sI l(SIR(x) > Pi) 



E 



J2 exp(-sJ' + PlhzWxW-")! 



Pit* 



> Pi 



(50) 



E e- s/ 'E ft 



E 



(51) 



E P [exp(-P(s(l + Pi) + P l Pr 1 M\ a ))} 



l + sPA\x\\ 



(52) 

where (a) follows from the fact that fading is independent 
of all the other random variables and (b) follows from the 
fact that h x ~ exp(l). Now, using the Laplace transform of 
I' and recalling r) = C{a) J2iLi MPlP'i > it can t> e further 
simplified to: 



E 



E 



eM-visil + PJ + PiP^Wxini) 



l + sB|b|| 



(53) 



and using Campbell Mecke theorem p3) to: 

/ exp(-,(,(14- ft H^-||, r) l) fc 

Jr2 l + sPi||x||- a 

With this we have now simplified both the terms of ( |48] l given 
respectively by ( |49] > and |54|. We now substitute the first term 
in (|46*1) and evaluate the integral with respect to s as: 



-l+2m/a 



exp ( -~r]s 2 / a ) ds 



?7 m a(m — 1)! 



, (55) 



where the solution follows from the substitution s 2 l a — > y 
followed by integration by parts. Now substituting the second 
term (given by |54l) in (|46l), we get the following integral: 



-l+2m/ a -r i {s(l+ji i )+0 i p- 1 \\x\\ o ') J i 



-dxdi 



1 + ,sPi||a;||- Q 

Now use the substitution (sPi) _1 / Q x — > x, which leads to 

" s -l+2m/a e — >7sa ((l+/3 i )+ l 8 i ||x|| 



(56) 



1 



(sPi) ° dxds. 

(57) 



Now exchange the integrals to obtain 



1 + \\x\\- a 



-dsdx. 
(58) 

Now the inner integral (with respect to s) can be evaluated 
directly using the definition of T(x) function or using the 
substitution s 2 l a s to obtain the below integral. 



XiPiPf ami 
2rfn-i 



dx 



(l + ||a;||-«)(l + A + AINII a )* (ro+1) 

(59) 



Now the above integral can be expressed as: 
1 f dx 



(1 + A)° (m+1) (l + ||a;|j-«)(l + T |^||x||«)S(" l + 1 ) 



Now using the substitution 1 + — , \\.r 



Pi 



(60) 

t~ l , the above 

expression can be simplified to 

2tt/37 2/q r(2m/a)r(l + 2/a) 
a(l + I3i) 2m / a T(l + (m + l)2/a) 
2 Fi(l, 2m/a, 1 + (m + l)2/a, (1 + ft)" 1 ), (61) 

where 2 Pi is the generalized hypergeometric function. Com- 
bining all the above we obtain the result. ■ 
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